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lHj; Abstract 

I ■ Fluids confined within narrow channels exhibit a variety of phases and phase transitions asso- 

^0 . ciated with their reduced dimensionality. In this review paper, we illustrate the crossover from 

I quasi-one dimensional to higher effective dimensionality behavior of fluids adsorbed within different 

' carbon nanotubes geometries. In the single nanotube geometry, no phase transitions can occur at 

. finite temperature. Instead, we identify a crossover from a quasi-one dimensional to a two dimen- 

^ I sional behavior of the adsorbate. In bundles of nanotubes, phase transitions at finite temperature 

•i-H . 

^ ' arise from the transverse coupling of interactions between channels. 
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I. INTRODUCTION 



One of the most interesting topics within modern condensed matter physics is that of 
phenomena in reduced dimensionahty, resuhing from some degree of spatial locahzation of 
the particles comprising a system, [ll] For example, chemists, materials scientists and physi- 
cists have created and explored numerous physical systems in which atoms and molecules 
are confined within quasi-one dimensional (QID) environments. The variety of these sys- 
tems is remarkable, such as the peapod geometry, i.e., a line of buckyballs within a carbon 
nanotube,jl, 0] fluids within artificial materials created by templating 0, 0] and QID op- 
tical lattices created by laser fields. 0] Unfortunately, as far as we know, there exists no 
comprehensive review of this general problem, although many relevant subfields have been 
summarized, jl, 0, [l^, lH The present paper addresses a small subset of this exciting re- 
search field. Specifically, we consider problems involving fluids, both classical and quantum, 
confined within QID channels, the focus of our group's research during the last decade, [l^ 

Here, the term QID refers to a system in which particles move in an external potential 
field ^(r) which is either constant or slowly varying in one direction (^), while V{y) is 
strongly localizing in the two other (transverse) directions. In the case of quantum particles, 
for which the transverse spectrum is discrete, one expects that the corresponding degrees of 
freedom are frozen out at low temperature (T); transverse excitation does occur at higher T, 
as determined by the gaps in the transverse spectrum of states. This plausible expectation 
is borne out in some cases, but we shall see that there can be dramatic consequences of the 
transverse degrees of freedom in other cases, even at low T. 

One of the many exciting aspects of strictly ID physics is its susceptibility to weak pertur- 
bations. The reason for this behavior arises from the fact that (in all practical situations) [l3| 
no phase transition can exist in a purely ID system at any finite T, even though the ground 
state may exhibit symmetry-breaking order. Thus, there do exist transitions at T identically 
equal to zero. A familiar example is the ID Ising model, for which the correlation length 
and susceptibility both diverge as T approaches 0. 

Since the purely ID system has no finite T transition, what brings about more interesting 
behavior in the QID case? As discussed below, the difference can arise from considering a 
set of parallel ID systems which are weakly coupled. Alternatively, the behavior can hap- 
pen because the system is only ID insofar as the transverse dimensions are finite, unlike the 

length in the z direction, L, which achieves the thermodynamic limit but the transverse 

dimensions are large enough to have an observable effect (e.g. low energy gaps). A third, 
more surprising, origin of interesting phenomena is when the system is a collection of nonin- 
teracting ID systems, with "quenched" heterogeneity, which are coupled to a particle bath, 
so they possess a common chemical potential. 

The next section discusses the conceptually and computationally simplest case of a QID 
system: a low density, noninteracting gas within a single channel; then, the solution of the 
one-particle Schrodinger equation determines the physical behavior. Section 3 considers the 
case of many interacting particles within a single channel; such a problem is often used 
as a model for fluids within regular or irregular porous materials. Section 4 considers the 
problem of QID channels containing fluids that interact with one another as well as with 
fluids occupying other channels. 
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FIG. 1: Left: Experimental heat capacity for ^He at various densities jnside ^He-coated FSM-16, 

The tube represents the 



from Ref. [5|]. Right: schematic depiction of the geometry, from Ref. 
environment experienced by a '^He atom. 



II. LOW DENSITY GAS IN A SINGLE QID CHANNEL 

Our first problem is conceptually and calculationally simple: a low density gas confined 
within a QID geometry; nevertheless, it provides interesting, and sometimes surprising, 
results. One example that has been studied extensively is that of a gas inside single carbon 
nanotubes, e.g. quantum fluids or the "peapod" case of Ceo molecules. [sl] Another is the case 
of quantum gases inside templated regular pores. ^ A third example is the so-called groove 
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region between two nanotubes, e.g. on the outside of a bundle of nanotubes. 

A classical noninteracting gas has a kinetic energy per particle of (3/2) ksT and a mean 
potential energy (U) determined by its interaction with the environment. For the case of 
a particle localized near the z axis within a channel, (U) = ksT, due to two transverse 
directions of excitation. Hence, the classical specific heat per particle is [C(T)/N]ciassicai = 
(5/2) /c^. For a quantum gas, instead, the transverse degrees of freedom are frozen out at 
low T, so one expects [C (T) / N]quantum = (1/2)^^. The generalization to D "effective" 
dimensions yields this expression for the dimensionless specific heat, C*, of a noninteracting 
Boltzmann gas: 

C* = C/iNks) = D/2 (1) 

Fig. [T] shows experimental results for C* in the case of ^He inside of FSM-16, a material 
consisting of straight hexagonal pores of cross-sectional distance of order 2 to 3 nm, pre- 
coated with a thin film of ^He. At the low densities shown here, the ^He gas can be assumed 
to be noninteracting, although interaction effects appear at higher density. [l3l The behavior 
observed inside FSM-16 can be understood by analogy to calculations in Fig. [2l for a 
noninteracting gas of ^He within a carbon nanotube. Note the overall similarity of these two 
figures. 

In Fig. [2],[l8| one observes that the low T limit is C* = 1/2, corresponding to a ID 
classical gas, as expected. At high T, instead, the limit is C* = 1; this limit is interpreted 
as that of a 2D gas moving on the inner surface of the nanotube. The bump at intermediate 
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FIG. 2: Heat capacity for a low density gas of ^He atoms within a single wall carbon nanotube, 
for cases R=0.8 nm (dashes) and 0.5 nm (full curve). Quantum statistical effects are omitted from 



the calculation. Adapted from Ref. 
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T is a general property found for spectra, such as that of the rigid rotor, 19| for which 
the inter-level spacing increases with quantum number. In the present case, the relevant 
spectrum is that arising from the azimuthal kinetic energy. 



Here z/ = 0, ±1, ±2, ... is the azimuthal quantum number and p is the radial coordinate, 
while is defined as a temperature characteristic of azimuthal excitation. The peak in the 
specific heat occurs at a temperature near 3 G for both radii considered in Fig. [2], so its 
position serves as a benchmark from which one can determine the value of (p). Note that 
inside a nanotube, the relevant value of (p) is typically R — a, where a is the gas-surface 
hard core interaction length. The key qualitative difference between the behaviors seen in 
Figs. [Hand [2] is that the experimental data in Fig. [1] plunge to C* = as T — 0. This is an 
effect of quantum degeneracy, manifested in Nernst's law, as found in explicit calculations 
which revise Fig. [2] by taking quantum statistics into account. j20| 

There have been many theoretical studies of gas adsorption in the presence of 
nanotubes. jl, 0, [13, 11, 12| In most treatments of these systems, one assumes that neigh- 
boring tubes are parallel. In that case, there exists a region of space- the so-called "groove"- 
which is a ID channel with a strongly attractive potential, created by the adjacent tubes. 
The adsorbed gas then exhibits QID behavior at low A^. However, one can also inquire 
about the case when the tubes are not quite parallel, but instead diverge, leaving a par- 
ticularly attractive region between them, with a minimum potential energy (Vq) located at 
equilibrium position tq. In a forthcoming study, [l|] we will report remarkable results for this 
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geometry, as exemplified in Fig. [31 The low T behavior is that of a gas in a ID harmonic 
potential, ^ Vq + kz{z — zoY/2, where kz is the force constant for particle motion parallel 
to the z axis, midway between and nearly parallel to the tubes' axes. We introduce the 
characteristic temperature for this motion, = h{kz/'my^'^ /ks and a reduced temperature 
T* = T/Tz- For T* << 1, the specific heat is of Arrhenius form, as seen in the figure, 
while for T* ~ 1, C* ~ 1, the specific heat of a 2D gas. That behavior might not have 
been anticipated, at first glance, because for small divergence half-angle 7, one might have 
expected ID behavior, i.e., C* = 1/2. Another surprise is the high T* limiting behavior, 
T* 7/4. This peculiar result arises because the relevant particles' motions are those in the 
plane perpendicular to the x — z plane of the nanotubes, for which the potential variation 
is unusual- proportional to y^. 

2.5r ■ ' ■ ■ n 



2.0- 




FIG. 3: Reduced specific heat C* for the gases He (dashed-blue) , H2 (solid-red) and Ne (short- 
dashed-black) between two nearly parallel nanotubes, each with R = 0.7 nm, and a divergence 
half-angle 7 = 0.5 degrees. From Ref. 



III. MANY INTERACTING PARTICLES IN A CHANNEL OF FINITE WIDTH 

From the perspective of phase transitions, a channel of finite width is a QID system, 
so that no true thermodynamic singularities can occur. This means that many attempts 
to explore behavior in porous media with single channel models (such as cylindrical and 
slit pores) cannot accurately describe phase transition behavior that is seen in genuinely 
3D porous media. Nevertheless, these models may provide good semiqualitative predictive 
power, sufficient for most purposes since the fully 3D geometry is not known. We have 
explored a variety of such models, with several different goals. These include assessing the 
accuracy of simplifying models, such as mean-field-theory (MFT) and the use of periodic 
boundary conditions. Both of these approximations are suspect, at first glance, due to the 
important role of fluctuations in ID statistical physics. 
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FIG. 4: Mean field isotherms compared with exact Monte Carlo isotherms as a function of reduced 
chemical potential and temperature, from Ref. 2l[. 



Fig. m presents adsorption isotherms for uptake in a single cylindrical pore, described by 
a lattice gas model. [2l[ In this model, continuous space is discretized into adsorption sites, 
which may be either occupied or empty in any microstate of the system. In the specific 
model used here, the set of sites (within one transverse section) consists of one axial site 
and seven "cylindrical shell" sites, corresponding to positions near the inner boundary of 
the nanotube. The ensemble of sites include an infinite sequence of such layers of eight sites. 
The energy of the system includes interactions between particles occupying these sites, plus 
interactions between particles and the substrate host. 

Fig. m compares (numerically) exact results with those obtained from MFT. The results 
are quite similar, overall, given the highly expanded scale of reduced chemical potential (/i*). 
Note that the spurious transition seen in MFT (at reduced temperature T = 0.5, /x* = —22) 
is not very different from the nearly discontinuous isotherm seen in the exact results. 



IV. REAL TRANSITIONS OF GASES WITHIN WEAKLY COUPLED QID 
CHANNELS 

Consider a geometry consisting of a set of parallel QID fluids, as in gases within or be- 
tween nanotubes comprising a bundle of such tubes. While no transition can occur for an 
isolated QID system, once coupling between such systems is present, a finite temperature 
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FIG. 5: Critical temperature as a function of the transverse interaction. Perturbation theory 
results and predictions from the anisotropic lattice gas model are compared. From Ref. Q]. 



transition can occur. Fig. [5] presents results for such a geometry. [3| The critical temper- 
ature Tc is shown as a function of the transverse coupling a, for the case of buckyballs 
confined within parallel nanotubes. As indicated, rather similar results were obtained from 
exact solutions for a lattice-gas model and from a perturbation theory. In the latter case, the 
unperturbed equation of state was the exact ID result for a model with Lennard- Jones inter- 
actions between the buckyballs and the perturbation was the weak van der Waals interaction 
between balls in adjacent tubes. 

Particularly striking in this figure is the singular behavior for small a, as blown-up in the 
inset. This behavior is well-known for the lattice-gas, for which the transition temperature 
satisfies (for weak coupling) ^2] 

^""^^ " /n(l/c) -/n[/n(l/c)] 

Here J/ is the longitudinal interaction and c = Jij Ji'vs, the anisotropy ratio of transverse to 
longitudinal interaction strengths. The singular behavior reflects the divergent susceptibility 
of the ID system at low T; as the correlation length diverges, larger regions of fluid in 
adjacent channels are coupled, so increases rapidly with increasing J^.The other notable 
result in this formula and the figure is that the characteristic energy and Tc scale are given 
by the longitudinal coupling, which is about 500 K for the peapod case, whereas naively one 
might have expected to be proportional to Jt. 
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